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Abstract. Quadrature formulas for f(x) dx where derivative terms need 
only be evaluated at a and b in the composite rule are identified. Error bounds 
are given when /: [a, b] — > R satisfies /( n_1 ) is absolutely continuous so that 
£ L p ([a,b]), and when /( n_1 ) is merely continuous. 



1. Introduction 

We are interested in finding numerical integration schemes for j f(x) dx where 
derivative terms in the composite rule need only be evaluated at the endpoints, 
a and b. Let p be a polynomial of degree n with leading coefficient 1/n! and 
/: [a, 6] — > R such that f( n_1 ) is absolutely continuous. Repeated integration by 
parts shows that 

f(x)dx = (b - (1^-^(1) f(b)-(b- 0,^-^(0)/ (a) 

(1.1) -{b - a)V"- 2 >(l)/'(&) + (b - a)V"- 2) (0)/'(a) 

+ ...-(-l)«(6-a)>(l)/(«- 1 )(6) 
+(-l)"(6 ~ a)>(0)/("- 1 )(a) + £ (o , 6) (p, /) 

where the error term is 

*W)(p, /) = (~l) n (6 - f P{x)f {n) {a + {b- a)x) dx. 

Jo 

This method is used in pQ. The composite rule is given in equation (3.1). Taking 
p(x) — x(x — l)/2 in equation (1.1) gives the trapezoidal rule. The Euler-Maclaurin 
formula follows from taking p(x) = B n (x)/n\ or p(x) — (B n (x) — B n )/n\ where 
B n (x) are Bernoulli polynomials and B n are Bernoulli numbers. The midpoint rule 
and Simpson's rule also follow from equation (1.1) [8]. 

We find the integration schemes where the derivative terms need only be evalu- 
ated at the endpoints of the interval [a, b] in the composite rule to be similar to the 
Euler-Maclaurin formula (Corollary 3.2). The most important of these integration 
schemes are essentially the Euler-Maclaurin formula. Similar formulas are derived 
in 0. 
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Section 2 outlines necessary facts on Bernoulli polynomials. The desired formulas 
are found in Section 3. In section 4 we give sharp error bounds when /(") e W{[a, b\) 
for 1 < p < oo and find asymptotic estimates for these bounds as n — > oo. In 
section 5 we weaken the assumption that is absolutely continuous to merely 

assuming that /t"" 1 ) is continuous. Then exists as a distribution and the 
error term can be evaluated with the distributional Denjoy integral. The error is 
estimated in terms of the Alexiewicz norm. 



2. Bernoulli Polynomials 
The Bernoulli polynomials can be defined recursively by B (x) = 1 and 

(2.1) -^B n {x) = nBn-^x) 
such that 

(2.2) \ B n (x)dx = 

Jo 

for n > 1. The Bernoulli numbers are defined as 

B n = B n (0). 

It then follows that 

(2.3) B n (x)=J2(^jB k x n - k . 

Integrating both sides of equation (2.3) on the interval [0, 1] leads to the equation 

£(T)*-o 

when n > 1. The Bernoulli numbers may then be calculated recursively with 
equation (2.4). 

It follows from induction that 

(2.5) B n (l-x) = (-l) n B„(a:). 

We may then conclude from equations (2.1), (2.2), and (2.5) that B n = when 
n > 1 is odd. It also follows from induction that B n (x) is monotonic on [0, 1/2] 
when n is even, and B n (x) does not change sign on [0, 1/2] when n is odd. 
Other properties of Bernoulli polynomials include 

(2.6) \B 2n (x)\ < \B 2n \ forze [0,1] 

(2.7) \B 2n (x)-B 2n \ < (2 - 2 1 - 2n )\B 2n \ for x e [0,1] 



1 - 2" +1 B n+1 
2™ n + 1 



,1/2 

(2.8) / B n (x)dx 
Jo 

(2.9) / B n (x)B m (x)dx = y > ' ' B m+n form,n> 1 
Jo {m + n)\ 

(2.10) ( _ 1) n-i ( 2 ^ 2 "^nW ^ CO s(27rx) as n ^ oo 

(2.11) ( 27r ) 2 " +lj32 "+ 1 ( x ) ^ sin(2 ^) as n -> oo. 
v ; v ' 2(2n+l)! v ' 
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The convergence in equations (2.10) and (2.11) is uniform on any compact set. All 
these results may be found in [3]. 

3. Main Theorem 

It is useful to define some new terms before proceeding further. Let p be a 
polynomial of degree n with leading coefficient 1/nl and / : [a,b] — > R such that 
y(n-i) j s aDSO l u tely continuous. Then 

I(a,b)(pJ) ■■= (b-afr^HVm-ib-a^-UiOMa) 

-(b - a)V"- 2 )(l)f (6) + (b- a)V"- 2 >(0)/'(a) 

+ ...-(-l)»(6-a)>(l)/("- 1 )(6) 

+ (-l)"(6-a)>(0)/("- 1 )(a) 

is said to be the integration scheme generated by p. We may now phrase equation 
(1.1) as 

rb 

f(x) dx = J( 0ib ) (p, /) + E (aM (p, /). 
Let h = (6 — a)/N for natural number N. The composite rule is 

(3-1) / /(*) ^ = Ifc tb) (p, /) + /) 

J a 

where 

JV-l 

I (l : b)(P^f) = X! / (a+'='!^+(fc+l)ft)(P, /) 
fe=0 

and 

N-l 

E (l,b)(Pif) = X! E (a.+kh,a+(k+l)h)(P> /)• 
fc=0 

In general, the derivative terms in the composite rule need only be evaluated at 
the endpoints, a and &, if and only if p^(0) = pW(l) for < £ < n — 2. In other 
words, this happens when the derivative terms of / form a telescoping series in 
the composite rule. This motivates us to say that p is a telescoping polynomial for 
I(a.b)(.Pi /) when p^ifi) = p( e \l) for < I < n — 2. In general, the function values 
of / in the composite rule cannot form a telescoping series since is linear. In 

the following theorem we charaterise the telescoping polynomials. 

Theorem 3.1. The polynomials 

BJx) . 



where c is an arbitrary constant, form the set of degree n telescoping polynomials. 
Proof. Let p be a polynomial of degree n with leading coefficient 1/nl. Write 



fc=0 v ' 



a k x „-k 



where ao = 1. Then 



w ^ fn - £ - k)W. 

k=0 y ' 
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Equating p (0) and p (1), it follows that p is a telescoping polynomial if and only 
if 

n-t-l 

y ^ = 

f-^ (n-£-k)\k\ 



for < £ < n — 2, i.e., 1 < n — t— 1 < n — 1. The a& are then defined recursively 
as in equation (2.4) for < k < n — 1 and a n is arbitrary. □ 

The composite rules produced by these polynomials are found by evaluating 
equation (3.1) with p(x) = B n {x)/n\ + c. 

Corollary 3.2. Let f : [a,b] — > R such that f( n ~ 1 '> is absolutely continuous. Then 
for each constant c, 



f 



f(x)dx = ^ w J(/(a ) + /(6 ) ) + ^^ r J^/(a + fc/ l ) 

k 



(3-2) +Et(^) (/^(a)-/^" 1 ^^ 

fe=2 ' ^ ' 

+c(fe - aTif^Ha) + f^Hb)) + E^ b _ a) (B n (x)/n\ + c, f). 

The fact that c is arbitrary begs the question of how c should be chosen. We 
will focus on the cases when c = —B n /n\ (n > 2) and when c = 0. The 
terms in equation (3.2) are not evaluated when c = —B n /nl. In this way, taking 
c = —B n /n\ gives the integration scheme with the least number of derivative terms 
generated by a degree n polynomial. The integration scheme generated by B n (x)/n\ 
is also generated by a polynomial of higher degree. This is not true when c ^ 0. 
Define p n {x) = B n (x)/n\ and q n (x) = B n (x)/n\ — B n /n\. These polynomials relate 
to the choices of c indicated above. When n > 3 is odd, p n = q n . Moreover, the 
integration schemes generated by p2n, P2n+i and q2 n +2 arc equal, and are essentially 
the Euler-Maclaurin formula. The following proposition shows some cases of when 
p n minimises the error estimate in equation (4.1). 

Proposition 3.3. (a) Choosing c = minimises \\p n + c|| 2 for each n > 0. 
(b) Moreover, c~0 minimises \\p n + c\\ r for each odd n, 1 < r < oo. 

Proof, (a) It follows from equation (2.2) that 

\\Pn + 4 2 2 = f [ Pn (x)} 2 dx + C 2 . 

Jo 

The desired conclusion is then easily observed. 

(b) The second assertion clearly follows from equation (2.5) when r = oo. When 

r < oo, 

,1/2 



f 1/2 r f 1 

\\p n + c\\ r r = / \p n {x) + c\ T dx + \p n (x)+c\ T dx 

JO Jl/2 

= / \Pn{x) + c\ r dx + / \p n (x) - c\ r dx 
Jo Jo 
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due to the asymmetry of p n about 1/2 when n is odd. The minimum of 
\pn{x) + c\ r + \p n {x) — c\ r occurs at c — 0. Hence, 

,1/2 

||rf=2 / \Pn(x)\ r dx<\\p n + c\\ r r . 



The inequality becomes strict for nonzero c when r > 1. □ 



4. L p Error Estimates 

We may use Holder's inequality to estimate the error terms under the assumed 
conditions. Let l/V + 1/.s = 1 for positive r, s. Choose / so that is absolutely 

continuous and G L s ([a,b]). Then 

(4-1) |^ (o , 6) (p,/)|<(6-ar+ 1 /'-||p|| r ||/(«)|| s 

where ||p|| r = (/„ x \p{x)\ r dx) ^ if r < oo, ||/W|| a = (f" \f^(x)\ s dx) 1 '* if s < oo, 

llplloo = max xe[0tl] \p(x)\, ||/ (n) ||oo = sup xe[aJl] \f^(x)\. Proposition 4.1 shows 
that this bound is sharp. 

Holder's inequality for series shows that 

N-l 

E ll/ (n) X[a+«»,-+(iH-l)/.]ll.<JV 1/r |l/ (n) ll. 
k=0 

when h = (b — a)/N. We may then estimate the error in the composite rules as 

(h — n\ n+1 / r 

(4-2) |*W)(P>/)| ^ HPllr||/ (n) || S . 

Proposition 4.1. The bound given in equation (4-1) is sharp. 

Proof. When r = 1, equality holds when / satisfies f( n \a+(b— a)x) — sgn(p(x)) for 
x G [0, 1]. When 1 < r < oo, choose / so that f( n \a+{b—a)x) = \p{x)\ r >s ; sgn(p(x)) . 
Then equality holds in equation (4.1) [5J Theorem 2.3]. Let {ipk} be the delta 
sequence given by tpk(x) = 1/fc for \x\ < l/(2fc) and V-'fc(a ; ) = otherwise, and 
let a; G [0, 1] maximise \p\ over [0, 1]. Choose fk to satisfy f^ (a + (6 — a)x) = 
V>fc(x — xo) if xo G (0, 1) and (a + (b — a)x) = 2ipk(x — xo) otherwise. Then 
ll/fc 111 = i b - a ) an£ l foP( x )fk n \ x ) dx -> p(x ) as fc -» oo. □ 

In general, calculating ||p|| r is a difficult problem. The following proposition 
gives some cases of when we are able to calculate this exactly for p n and q n . 

Proposition 4.2. (a) \\q n \\i = \B n \/n\ when n is even. 

(b) \\Pn\\i = (1 - 2" +1 )|g„+i|/((n + 1)2"- 1 ) when n is odd. 

(c) \\p n \\ 2 = v /|B 2n |/(2n)r and [| 9n [[ 2 = ^\B 2n \/{2n)\ + {B n /n\y '. 

(d) Whennis even, \\p n \\ x = \B n \/n\ and Wq^ = |g„(l/2)| < (2 - 2 1 - n )\B n \/n\. 

Proof, (a) Let n be even. Then q n does not change sign on [0, 1] by equation (2.6). 
Hence, 



\q n (x)\dx 



l 



q n (x) dx 



o 



= \B n 
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(b) Let n be odd. Then p n is asymmetric about 1/2 while not changing sign on 
[0,1/2]. Hence, 

ri/2 

q n (x) dx 



J 

Jo 



\p n (x)\ dx = 2 



l-2" +1 |B n+1 | 



n+ 1 



The final equality follows from equation (2.8). 

(c) The first inequality is evident from equation (2.9). The second then follows 
from the first while noting equation (2.2). 

(d) The first equality is essentially equation (2.6). The second equality comes 
from q n being monotonic over [0, 1/2] while being symmetric about 1/2 with q n (0) = 
0. The inequality is essentially equation (2.7). □ 

We now shift our attention to finding asymptotic estimates of \\p n \\r and ||g n || r 
as n — > oo. Since the convergence in equations (2.10) and (2.11) is uniform, 

r r i 



(4.3) / \p 2n {x)\ r dx 
Jo 

(4.4) f \p 2n+1 (x)\ r dx 
Jo 

(4.5) / \q 2n (x)\ r dx 
Jo 



(2 



2 V f L 

,„ ) / I cos(27nc)| r dx as n — > oo 
*? n ) Jo 



2n+l 



(27T) 



(2tt) 2 " 



/ | sin(2 
Jo 



TTx)\ r dx as n — y oo 



cos(27rx) + l| r dx as n — > oo 



and 
(4.6) 

(4.7) 



IK I 



l<?2n| 



(2tt)' 
4 



as ji — t oo 



as n — > oo. 



(2tt) 2 " 

Evaluating equations (4.3), (4.4) and (4.5) at r = 1 gives: 



(4.8) 
(4.9) 

It is well known that 

r n/2 

(4.10) 



IK Hi 

ll<72n||l 

; dx - 



(277)"+! 
2 



as 7i — y oo 



(27T) 



2d 



as 7i — y oo. 



" /2 r J l-3-5---(r-l)7r 

cos xdx — — - — ■ — — 

2-4-6--T 2 



/ /2 r , 2-4-6---(r-l) 
sm xdx — / cos xdx = 



when r > 2 is even and 

f 7r/2 

(4.11) , ... , 

1 1 1-3-5--T 

when r > 3 is odd. Using equations (4.10) and (4.11) we find the following asymp- 
totic estimates for \\p n \\r and \\q n \\r when r > 2 is an integer: 



pit j Z pT. 

/ sin r xdx = 
Jo Jo 



(4.12) 



\\Pn\ 



(2tt)™ 



l-3-5---(r-l) 
2-4-6--T 



l/r 



as n — >■ oo 



when r > 2 is even, and 



(4.13) 



IK I 



2r+l 



(27r) r 



2-4-6---(r- 1) 2 
1-3-5--T ^ 



l/r 



as 7i — y oo 
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when r > 3 is odd. 

5. Generalised Error Bounds 

We initially assumed that /( n_1 ) was absolutely continuous on [a, b]. We may 
weaken this assumption to merely assuming /( n— x ) is continuous if we evaluate 
integrals in the distributional Denjoy sense. The distributional Denjoy integral 
allows integrating derivatives of any continuous functions where derivatives are 
taken in the distributional sense. 

Define the set of test functions to be T>([a, b]) — {(f) E C°°([a,b]) \ 4> has support in (a,b)} 
where the support of a function is the closure of the set of points where the function 
is nonzero. A sequence {</>„} C V([a,b]) is said to converge to 4> s V([a,b]) when 
there exists a compact set K C (a, b) so that </>„ has support in K for each n and 
4>n — ► 4>^ m ^ uniformly. The distributions are the continuous linear functionals 
on the space of test functions; the set of distributions is denoted V([a,b]). Let 
/ : [a,b] — > R be an integrable function. Then (f,<f>) — J f(x)(f>(x) dx defines a 
distribution. 

Let T be a distribution. Then T so that (T' , 4>) = — (T, 4>') for each 4> 6 T>([a, b\) 
is said to be the distributional derivative of T. If T is taken to be an integrable 
function / : [a, b] — > R with a pointwise derivative at xq then that pointwise value 
may be identified by evaluating distribution / with a delta sequence of test functions 
about xq- 

The definition of the distributional Denjoy integral is simple. If / is the distri- 
butional derivative of a continuous function F : [a, b] — > R then J / = F(b) — F(a) 
defines the integral over a compact interval. Take / to be the distributional deriv- 
ative of a continuous function F : [a,b] —> R and let g : [a,b] —> R be of bounded 
variation. The integration by parts formula for the distributional Denjoy integral 
is given by 

(5.1) ( fg = F(b)g(b)-F(a)g(a) - [ Fdg. 

In this way, equations (1.1) and (3.1) still hold with this weakened assumption 
that /t'^ 1 ) be merely continuous. The Beesack-Darst-Pollard inequality for the 
distributional Denjoy integral is given by 



(5.2) 









f fg 


< 




J a 




J a 



inf \g\ + ||/|U e V[ 0) 5]fl. 

[a,b] 



Here || ■ \\a c is the Alexiewicz norm given by 

||/|U= max \F(x)-F(y)\ 

and V[ a Mg is the total variation of g. All this on the distributional Denjoy integral 
is found in [B]. 

Applying equation (5.2) to E/ a ^\(p n + c, f) gives the estimate 
(5-3) \E ia , b) (p n + cJ)\ < (b-a) n \\fW\UV [0 , 1]Pn 

when p n + c has a real root in [0, 1] . This is true for both p n and q n . It then follows 
that 

(5-4) \E^ a) (Pn + c,f)\ < { ±^\\fM\\ Ac V [0A]Pn 
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when p n + c has a real root in [0, 1]. 

The total variation of p n may be calculated from the maxima and minima of 
Bernoulli polynomials on [0,1]. See [1] for estimates of these upper and lower 
bounds. 

We use the uniform convergence in equations (2.10) and (2.11) to give the as- 
ymptotic estimate 

g 

v [o,i\Pn ~ j^y^ as n -> oo. 
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